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Abstract 

An extended quark mass density- and temperature- dependent model which includes the cou- 
plings between quarks and the cr-mesons, w-mesons is suggested. The MIT bag boundary constrain 
has been given up and the interactions between quarks and mesons are extended to the whole free 
space. We show that the present model is successful to describe both the saturation properties and 
the deconfinement phase transition of nuclear matter. When the effective nucleon masses vanish 
and the bag radius tends to infinite, the quark deconfinement phase transition takes place. The 
corresponding QGP phase diagram is addressed. 
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I. INTRODUCTION 



Because of the non-perturbative properties of QCD at low energy regions, it is very 
difficult to study nuclear system by using QCD as a fundamental theory directly. Various 
phenomenological models based either on hadron degree of freedom or quark degree of 
freedom have been suggested. The quark-meson coupling (QMC) model, first proposed by 
Guichon[l|, which describes nuclear matter as a collection of non-overlapping MIT bags, 
scalar a meson and vector u meson is one of such successful candidates. The quarks inside 
the MIT bag couple with the scalar a meson and vector uj meson self-consistently. By means 
of this model and the mean field approximation, many dynamical and thermal properties 
of nucleon systems and hyperon systems have been studied both at zero temperature and 
finite temperature^, 3, fl U fl Q, 3 

Although QMC model is successful for describing the physical properties of nuclear sys- 
tem, many shortcomings arise when one use this model to discuss the quark deconfinement. 
The first difficulty comes from that this is a permanent quark confinement model. The 
nucleon corresponds to a MIT bag and the mechanism of quark confinement in MIT bag 
model being that the normal flow of quark current is zero at the bag surface. This boundary 
condition can not be changed by temperature and density. The second difficulty arises from 
the many-body calculations. If we hope to do the nuclear many-body calculations beyond 
mean field approximation by quantum field theory, it is essential to find the propagators of 
quark, a meson and u meson respectively. But the constrain of MIT bag boundary condi- 
tion and the interaction between quarks and mesons limited within the bag regions present 
obstacles to get the corresponding propagators in free space. 

On the other hand, another effective model, namely, the quark mass density- depedent 
(QMDD) model, was suggested by Fowler, Raha and Weiner tj. According to this model, 
the masses of u, d quarks and strange s quark (togather with the corresponding anti-quarks) 
are given by 

m q = — (q = u,d,u,d) (1) 
n Q 

B 

m s ,s = m s0 H (2) 

where B is the vacuum energy density inside the bag and uq is the quark number density. 
The basic hypothesis Eqs.(l) and (2) corresponds to a quark confinement mechanism because 



if quark goes to infinite, the volume of the system tends to infinite, tiq approaches zero 
and m„ approaches to infinite. The infinite quark mass prevents the quark from going to 



r 

infinite jlOl This confinement mechanism is very similar to that of the MIT bag model. 
As was shown in ref. Q], the properties of strange matter in QMDD model are nearly the 
same as those obtained by the MIT bag model. Although the confinement mechanism is 
similar, the advantage of QMDD model is that it does not need to introduce a boundary 
condition to confine quark as that of the MIT bag model. 

As was pointed by refs. [12I . k| , if we use QMDD model to investigate the thermodynam- 
ical properties of nuclear system at finite temperature, many difficulties will emerge. For 
example, it cannot reproduce a correct lattice QCD phase diagram qualitatively because 



the masses of quarks, and then the temperature tends to infinite when — > 



121. The 



reason is that the confinement mechanism of QMDD mode 



is still permanent. To over- 
come this difficulty, in a series of our pervious papers we suggested a 
quark mass density- and temperature- dependent (QMDTD) model. Instead of the perma- 
nent MIT quark confinement mechanism, we employed the Friedberg-Lee (FL) soliton bag 
model. Since the confinement mechanism of FL model comes from the interaction between 
quarks and a nontopological scalar soliton field, and the spontaneously broken symmetry 
of the scalar field will be restored at a finite temperature, the soliton solution will disap- 
pear and the quark will deconfine at the critical temperature. In this nonpermanent quark 
confinement model, the vacuum energy density B equals the different value between the per- 
turbative vacuum and physical vacuum, and it depends on temperature. Instead of Eqs.(l) 
and (2), we introduced 
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in QMDTD model. The masses of quarks not only depend on density but also on temper- 
ature. Obviously, the QMDTD model reduces to QMDD model at zero temperature. By 
means of the QMDTD model, the physical properties and the stability of strangelets|13], 
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181 ] at finite temperature have been studied 



the dibaryon system |15( and the strange star 
and the results are fine. 

Although the masses of quarks reflect the confinement characteristic in the QMDTD 
model, it is still an ideal quark gas model. If we hope to investigate the saturation properties 
and the deconfinement phase transition of nuclear matter, the quark-quark interaction must 
be taken into consideration. In ref.Q], a quark and non-linear scalar field coupling is 
introduced to improved QMDD model at zero temperature. We proved that many physical 
properties given by FL soliton bag model can be mimicked by the improved QMDD model. 

This paper evolves from an attempt to extend our study to finite temperature. We will 
introduce the scalar a meson, vector uj meson and the couplings between quarks (u, d) 
and u meson, uo meson in the extended quark mass density- and temperature- dependent 
(EQMDTD) model. Our EQMDTD model is similar to that of Walecka model and QMC 
model. The basic differences between the EQMDTD model and the Walecka model are: 
(l)we replace the nucleon in Walecka model by quark in EQMDTD model; (2)rather than the 
structureless point-like nucleon in Walecka model, the nucleon corresponds to a "cluster", 
say, bag which consists three quarks with temperature- and density- dependent masses 
given by Eqs(5) and (6). Besides, the differences between the EQMDTD model and the 
QMC model are: (l)instead of the MIT bag in QMC model, we incorporate quarks with 
density- and temperature- dependent masses in EQMDTD model; (2)in place of the quark-u 
quark-o; interactions restricted within the bag region in QMC model, the quark-cx, quark- 
uj interactions are spreaded to the whole free space, because the constraint of MIT bag 
boundary condition is discarded. Moreover, since our model are based on the FL model, the 
quark deconfinement phase transition can take place, we will address the behavior of the 
saturation and the deconfinement phase transition for EQMDTD model in this paper. 

The organization of this paper is as follows. In the next section, we will give a brief 
description of EQMDTD model and the main formulae. Our numerical results will be 
presented in Section III. The last section is devoted to discussion and conclusion. 
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II. THE EQMDTD MODEL 



The Lagrangian density of the EQMDTD model is 

£ = i> -m q + {g%a - 9^,)] ^ + ^d,ad»a - l -m 2 a a 2 - -F^ + ^m>^ (9) 

where the quark mass m q is given by Eqs.(5) and (6), m a , are the masses of a and 
uj mesons, = d,u u — dyU)^ g% and g% are the couplings between quark-a meson and 
quark-cu meson respectively. We neglect the s quark in the following discussion. 
The equation of motion for quark field is 

[y {id, + gtuo,) - (m, - gla)] ^ = (10) 

Under mean field approximation, the effective quark mass m* is given by 

m* q = m q - g q G a (11) 

In nuclear matter, three quarks constitute a bag, say nucleon, and the effective nucleon mass 
is obtained from the bag energy and reads 

M* N = E bag = Y.E q = |W £ £ £,(*)(/,(*) + /,(*)) (12) 

1 q=u,d k 



where e q (k) = y'm* 2 + k 2 is the single particle energy, 7 is the degeneracy, f q {k) and f q (k) 
are the Fermi distributions of quark and anti-quark respectively 

f*W = e ^ q {k)-^ q ) + 1 ( 13 ) 

= e P(e q (k)+„ q ) + I ( 14 ) 

The chemical potential ji q of quark is given by 

n Q = 7 E £(/#) - /,(*)) (15) 

q=u,d k 

3 = ^nR 3 n Q (16) 
The bag radius R is determined by the equilibrium condition for the nucleon bag 
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In nuclear matter, the total energy density is given by 

2 



£ 



matter 



E 



i=N,P ^) 

where ps is the density of nuclear matter 



<?ky/M? + k* (TH(k) + n % {k)) + -fL 2 s + \mla 2 

ZiTTl., / 
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(19) 



n»(A;) and 7ij(/c) are the Fermi distributions of nucleon and anti-nucleon respectively, 

rii(k) 
fii(k) 



1 



/ij is the chemical potential of the ith nucleon. In Eq.(jI8j). g u is the cou 
and to meson and it satisfies g u = 3g«. As that of QMC modelQ, S, 0, Q,H fl 3 
mean field of a meson is determined by self-consistency condition 

5E ma n er 



(20) 
(21) 
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which yields 



a 
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ing between nucleon 
the scalar 

(22) 



Mf + k 2 \ 9(7 



{riiik) + Wi(fc)) 



(23) 



Eqs.(5),(9)-(20) form a complete set of equations and we can solve them numerically. Our 
numerical results will be shown in the next section. 



III. RESULTS 

Before numerical calculation, let us discuss the effective quark mass m* carefully because 
this quantity affects our results directly. To take the medium effects into account, Jin and 
Jinnings suggested a modified QMC model [3] in which the bag parameter B depends on den- 
sity. In EQMDTD model, based on FL model, the bag parameter B depends on temperature 
(Eq.(3)), and the quark mass m q depends on quark density uq. As was shown in section I, 
this quark density dependence of m q corresponds to a quark confinement mechanism only. 
The medium effect has not yet been considered. To exhibit the medium effects, instead of 
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B = B(uq) in the modified QMC model, we introduce an ansatz that the effective coupling 
of quarks and a meson g q a is a function of quark density uq, which can be expanded as 

9l=gfVn Q + gl m lnl + gf)/nl (24) 

Eq.(9) becomes 

. B„ [l - (T/T c f] - (g?» + 9 »W/n g + gf>/,, 2 Q )a 
m , = (25) 

In fact, the density- dependent couplings of NNw, NNp : NNa have been employed by 
many authors to discuss many different problems Q, fl H Q • 

Eq. ()24jl is an extension 

to quark level only. The values of the adjusted parameters g^°\ and g q ^ will be 
determined below. 

Now we discuss the parameters in EQMDTD model. Firstly, we choose = 783 MeV, 



m„ = 509 MeV as that of refs. 
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25]; T c = 170 MeV as the deconfinement temperature at 
zero baryon density. To fix the nucleon mass = 939 MeV, we take B = 173 MeVfm -3 . 
Beside these parameters, there are still four parameters, namely, g%>°\ g a ^ and g q ^ 
needed to be fixed. Obviously, the behavior at the saturation point and the deconfinement 
phase transition must be explained if EQMDTD model is successful. We see below that the 
quark deconfinement phase transition will take place at the point — > and bag radius 
R — > oo. Therefore, we fix parameters g$°\ g^ 1 ' and g q ^ by the condition that at zero 
temperature, binding energy E = —15 MeV, the ratio of the effective nucleon mass to the 
nucleon mass M^/Mjv = 0.6 at the saturation density po — 0.15 fm~ 3 , and the deconfinment 
phase transition takes place at ps = 10po- We find 

gf> = 5.93/m" 3 , gf^ = -0.747/m" 6 , gf^ = 0.035/m" 9 , g* = 4.23 

Now we are in a position to do numerical calculations and the results are shwon in Fig. la- 
Fig. 9b. In Fig. la, we omit the contribution of cr-meson and tu-meson and depict the bag 
energy as a function of bag radius at different temperatures T = 0, 30, 60, 100, 150, 170 MeV 
respectively. In this specific case, the EQMDTD model reduces to the original QMDTD 
model without quark- meson interaction. Fig. lb is the same diagram as that of Fig. la but 
that contribution of meson fields has been taken into consideration. We fix a = 20 MeV. It is 
found from Fig. la and Fig. lb that the radius of the bag given by the minimum of the curve 
increases with increasing temperature. When T < lOOMeV, the bag radius varies merely 



slightly with increasing temperature. This result is similar to that of the QMC model or 
the FL model. However, when the temperature approaches to the critical temperature t c , 
the bag radius increases significantly. 

Noting that the a connects to the baryon density pb by Eq.(23) directly, instead of ps, 
we employ a to show our results first. The curves of effective nuclear mass M* vs. the mean 
field a at different temperatures T = 0, 30, 60, 100, 130, 160 MeV are shown in Fig.2. We 
see that the effective nucleon mass decreases when a increases. In particular, we find a very 
interesting phenomenon: when a approaches to a fixed point <r c , M* drops to zero rapidly. 
This fixed point a c decreases as the temperature increases. The result M* = means that 
the nucleon is dissolved and the quark deconfinement phase transition takes place. To show 
this phase transition more clearly, we draw the curves R vs. a at different temperatures 
T = 0, 30, 60, 100, 130, 160 MeV in Fig.3. We see that at the same critical point a c shown 
by Fig.2, the bag radius R tends to infinite rapidly. This result confirms that this is indeed 
a critical point of deconfinement phase transition, because at this case quarks occupy the 
whole system. Similarly, for fixed a = 0,20,40,60,70,75 MeV, we show M* vs. T curves 
and R vs. T curves in Fig.4 and Fig.5 respectively. We see that for every curve, M* drops to 
zero and R tends to infinite at a critical temperature t c . The critical temperature t c increases 
when a c decreases. There results exhibit that the quark deconfinement phase transition can 
be described by EQMDTD model naturally. 

Now let us show the connection between a c and the baryon density ps- Employing Eqs. 
(19) and (23), one can obtain the relation of a and ps- The result is shown in Fig. 6. For 
a fixed temperature, a is a monotonical function of pB- It means that a can play the same 
role as that of ps- For example, M* vs. o curve has the same behavior as that of M* vs. pb 
curve. The curve in Fig. 7 shows that when a c increases, the critical temperature t c decreases. 
It corresponds to the result that when p c B increases, t c decreases for the quark deconfinement 
phase transition. This is of course very reasonable. In fact, Fig.7 is the deconfinement phase 
diagram for EQMDTD model. It is very similar to that given by lattice calculation. 

Now we turn to address the saturation properties of nuclear matter at low temperature 
and low density by using the EQMDTD model. Our results are shown in Fig. 8, Fig.9a 
and Fig.9b. The curves of £ ma tterl Pb — M N vs. ps for different temperatures are shown in 
Fig. 8 where the curve at T = MeV refers to the saturation curve and the corresponding 
minimum is the saturation point. The changes of the effective nucleon mass M* to ps are 
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shown in Fig. 9a and Fig.9b where the curves in Fig. 9a limit to the low density region (0 
fm~ 3 < p B < 0.20 fm~ 3 ), and curves in Fig.9b cover a large density region. We find that 
M* decreases when density and/or temperature increases. This result is similar to the of 
QMC model[47] and Walecka model. At the saturation point, M* /M = 0.6, as indicated. 

IV. SUMMARY 

In summary, we present an extended quark mass density- and temperature- dependent 
model which is motivated by the QMDTD model and the quark meson coupling model. The 
a and uj fields are assumed to couple with the u, d quarks. The MIT boundary constrain 
has been given up and the non-permanent quark confinement mechanism is mimicked by 
the quark mass density- and temperature- dependence. It is shown that under mean field 
approximation this model provides a reasonable description of the saturation properties 
for nuclear matter in the region of low temperature and density. While at sufficient high 
temperature and/or density, the deconfinement phase transition takes place automatically. 
The phase transition is characterized by that the bag radius tends to infinite and the effective 
baryon mass M* vanishes simultaneously. The deconfinement phase diagram which has 
the same behavior as that of lattice QCD calculation is addressed. We emphasize that in 
our calculation the deconfinement process is treated consistently within one unified model. 
This is different from the usual investigation which employs two separate models to describe 
respectively the hadronic phase and quark gluon plasma phase, and use the Gibbs conditions 
to determine the phase diagram[26,27]. 
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Fig. 1a 



FIG. 1: The bag energy as a function of bag radius without meson at different temperatures 
T = 0, 30, 60, 100, 150 and 168 MeV, with a = MeV, the radius for nucleon is determined by the 
minimum of bag energy 1 1 
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Fig. 1b 



FIG. 2: The bag energy as a function of bag radius without meson at different temperatures 
T = 0, 30, 60, 100, 150 and 170 MeV, with a = 20 MeV. 
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FIG. 3: Effective nucleoli masses vs. the a at different temperatures T = 0, 30, 60, 100, 130 and 
160 MeV respectively i o 




FIG. 4: Nucleon radius vs. the a at different temperatures T = 0, 30, 60, 100, 130 and 160 MeV 
respectively i a 




Fig.4 



FIG. 5: Effective nucleoli masses vs. the temperatures at different values a = 0, 20, 40, 60, 70 and 
75 MeV respectively i cr 




FIG. 6: Nucleon radius vs. the temperatures at different values a = 
respectively i c 



0, 20, 40, 60, 70 and 75 MeV 




1 1 1 1 1 1 r 

0.0 0.5 1.0 1.5 

P B (fm 3 ) 



Fig.6 



FIG. 7: a vs. baryon density ps at different temperatures T = 0, 50, 100 and 130 MeV respectively 
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Fig.8 



FIG. 9: Saturation curve of nuclear matter at different temperatures T = 0, 100 and 130 MeV 
respectively i q 




Fig.9a 



FIG. 10: Effective nucleon mass vs. baryon density at different temperatures T = 0, 50 and 100 
MeV respectively 2D 
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Fig.9b 



FIG. 11: The same curve as Fig. 4b with a wider range of baryon density at different temperatures 
T = 0, 50 and 100 MeV respectively 2 1 



